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Abstract 

Computational  design  analysis  of  a  high  performance  cathode  is  a  cost-effective  means  of  exploring  new  microstructure  and  material 
options  for  solid  oxide  fuel  cells.  A  two-layered  porous  cathode  design  has  been  developed  that  includes  a  thinner  layer  with  smaller 
grain  diameters  at  the  cathode/electrolyte  interface  followed  by  a  relatively  thicker  outer  layer  with  larger  grains  at  the  electrode/oxidant 
interface.  Results  are  presented  for  the  determination  of  spatially  dependent  current  generation  distributions,  assessment  of  the  importance 
of  concentration  polarization,  and  sensitivity  to  measurable  microstructural  variables.  Estimates  of  the  electrode  performance  in  air  at 
700  °C  indicate  that  performance  approaching  3.1  A/cm2  at  0.078  V  is  theoretically  possible.  The  limitations  of  the  model  are  described, 
along  with  efforts  needed  to  verify  and  refine  the  predictions.  The  feasibility  of  fabricating  the  electrode  configuration  is  also  discussed. 
©  2003  Published  by  Elsevier  B.V. 
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1.  Introduction 

A  solid  oxide  fuel  cell  (SOFC)  is  composed  of  a  dense 
electrolyte  sandwiched  between  porous  electrodes.  The  key 
electrochemical  reactions  occur  mostly  on  the  surfaces  at  the 
electrolyte/electrode  interfaces,  thus  enabling  the  harvesting 
of  electrons  in  an  electrical  circuit  to  produce  useful  power. 
The  dense  electrolyte  conducts  oxygen  ions  from  the  cath¬ 
ode  to  the  anode,  but  prevents  mixing  of  the  fuel  and  oxidant 
in  the  gas  phase,  where  electrons  cannot  be  harvested.  The 
porous  electrodes  permit  passage  of  the  gases  to  the  elec¬ 
trolyte/electrode  interfaces  where  the  reactions  occur.  A  rel¬ 
atively  thick  porous  anode  is  often  used  to  provide  structural 
support  for  the  assembly.  The  cathode  is  often  relatively  thin 
to  minimize  its  contribution  to  the  overpotential  by  polar¬ 
ization  losses.  It  is  generally  recognized  that  a  significant 
portion  of  these  polarization  losses  originate  in  the  cathode. 
Consequently,  much  attention  has  recently  been  focused  on 
improving  cathode  performance  through  two  means:  im¬ 
proved  materials  and  improved  microstructural  designs.  Im¬ 
proved  materials  include  mixed  ionic  electronic  conductors 
(MIECs),  which  essentially  increase  the  surface  area  avail¬ 
able  for  conversion  of  gaseous  oxygen  molecules  to  oxygen 
ions.  Such  a  material  has  been  developed  in  our  laboratory 
[1],  and  was  employed  in  this  work.  The  present  paper  fo¬ 
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cuses  on  methods  to  improve  the  microstructural  design  of 
the  cathode. 

An  objective  of  this  work  is  to  design  a  solid  oxide 
fuel  cell  (SOFC)  cathode  exhibiting  a  low  area  specific 
resistance  (ASR  =  0.1  £2  cm2)  and  capable  of  high  current 
output.  Such  an  objective  can  be  attained  best  by  a  coupled 
experimental-modeling  approach:  the  modeling  helps  to 
guide  the  experiments  and  the  experiments  provide  data  for 
fitting  the  parameters  of  the  model.  Since  modeling  is  often 
cheaper  than  a  long  series  of  Edisonian  experiments,  costs 
are  often  minimized  with  this  approach.  This  paper  de¬ 
scribes  modeling  efforts  and  includes  selected  experimental 
data,  which  are  reported  in  more  detail  elsewhere  [1]. 

An  initial  step  was  to  review  the  existing  models  in  the 
literature,  and  select  an  approach  that  was  both  pragmatic 
and  thorough,  in  addition  to  exhibiting  direct  linkages  with 
the  experimental  data.  We  found  four  categories  of  models, 
which  are  described  below  in  terms  of  complexity  and  in¬ 
put  data  requirements.  The  following  paragraphs  are  not  an 
exhaustive  review,  and  contain  only  representative  examples 
of  each  model  type. 

In  the  first  category  are  detailed  models  by  Svensson 
et  al.  [2].  This  model  treats  the  classic  three  phase  bound¬ 
ary  (TPB)  problem  explicitly  by  addressing  the  individual 
mechanisms  involved  (surface  adsorption,  dissociation,  elec¬ 
tronic  exchange,  surface  diffusion).  Non-linear  second  order 
differential  equations  are  derived  and  solved  numerically, 
with  the  proper  boundary  conditions  (six  are  required).  Al¬ 
though  the  model  certainly  contains  enough  technical  depth, 
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it  also  requires  the  estimation  of  about  a  dozen  parameters 
(such  as  surface  diffusion  coefficients)  that  generally  have 
uncertainties  of  two  to  four  orders  of  magnitude.  Estima¬ 
tion  of  these  parameters  would  require  an  extensive  series  of 
well-controlled  experiments,  at  appreciable  cost.  This  model 
was  not  suitable  for  the  present  program.  Similar  conclu¬ 
sions  have  been  reached  elsewhere  [3]. 

In  the  second  category  are  models  based  on  a  homog¬ 
enized,  effective  medium  concept  for  the  porous  cathode 
material  [4,5].  The  TPB  mechanisms  described  above  are 
lumped  into  a  single  parameter  called  the  surface  exchange 
coefficient  (7TS,  cm/s),  which  is  experimentally  measurable 
on  a  routine  basis,  e.g.  [6].  A  significant  database  exists 
in  the  literature  for  a  range  of  mixed  ionic  electronic  con¬ 
ductors  (MIEC)  that  are  important  in  the  present  effort.  In 
this  respect,  the  models  in  [4,5]  are  suitable  because  of  the 
direct  connection  to  experimental  tasks.  However,  two  fac¬ 
tors  present  problems.  The  first  was  the  mathematical  com¬ 
plexity  resulting  from  treatment  of  four  polarization  terms 
(Butler- Volmer,  chemical  resistance  caused  by  the  exchange 
process,  concentration  polarization,  and  gas  capacitance  at 
the  cathode/gas  interface).  Second,  communication  with  the 
author  [4,5]  indicated  that  several  typographical  errors  ex¬ 
isted  in  the  literature,  thus  requiring  re-derivation  of  the 
model  to  ensure  correctness.  The  investment  required  for 
this  type  of  model  also  rendered  it  unsuitable  for  the  present 
work. 

In  the  third  category  are  models  by  Deng  et  al.  (DZA) 
[7]  (without  concentration  polarization),  and  [8]  (with  con¬ 
centration  polarization).  These  models  use  the  same  surface 
exchange  concept  as  in  [4,5],  but  treat  only  the  chemical 
resistance  term  [7]  and  the  concentration  polarization  term 
[7,8].  Numerical  input  requirements  are  comparatively  mod¬ 
est,  and  the  models  are  mathematically  tractable.  Although 
solutions  are  available  in  closed-form,  these  solutions  are 
for  a  homogeneous  material,  so  that  spatially  variant  mate¬ 
rial  properties  are  not  treated.  Consequently,  numerical  solu¬ 
tions  of  the  differential  equations  will  ultimately  be  required 
for  multilayered  cathode  concepts.  The  models  are  also  well 
suited  for  the  high  conductivity  materials  used  in  this  inves¬ 
tigation,  and  were  selected  for  the  calculations  herein. 

In  the  fourth  category  are  the  more  classical  models, 
e.g.  [9],  which  generally  treat  all  polarization  terms  in  a 
semi-empirical  manner.  They  require  input  of  a  modest 
amount  of  experimental  data,  most  of  which  is  available, 
and  are  mathematically  tractable.  However,  linearizations 
in  the  model  have  apparently  caused  the  loss  of  important 
effects,  particularly  for  gas  diffusion  through  the  porous 
material.  The  manifestation  of  this  problem  was  the  rather 
high  values  derived  for  the  tortuosity,  indicating  that  several 
unknowns  and  uncertainties  had  numerically  accumulated 
in  this  parameter.  This  type  of  model  was  thus  not  suitable 
for  the  present  effort. 

We  selected  the  models  by  Deng  et  al.  (DZA)  [7,8]  for 
this  work  because  of  their  direct  relationship  to  experimen¬ 
tal  measurables,  treatment  of  the  two  most  important  polar¬ 


ization  terms,  mathematical  tractability,  and  the  availability 
of  closed-form  solutions  for  initial  estimates.  The  following 
paragraphs  describe  preliminary  results  using  these  models. 

2.  Engineered  cathode  development 

2.7.  Estimation  of  current  distribution 

In  the  early  stages  of  this  investigation,  the  lanthanum 
strontium  copper  ferrite  (LSCuF)  cathode  material  being 
tested  exhibited  very  low  ASR  and  high  current.  It  was  there¬ 
fore  thought  that  concentration  polarization  was  negligible 
up  to  about  0. 1  V.  The  DZA  [7]  model  was  extended  to  per¬ 
mit  estimation  of  possible  oxygen  potential  gradients,  and 
used  to  estimate  the  current  distribution  throughout  the  cath¬ 
ode  material.  Model  parameters  are  defined  in  Appendix  A, 
and  the  model  details  are  given  in  Appendix  B,  along  with 
a  benchmark  case  for  LSCoF. 

Preliminary  results  for  the  very  efficient  mixed  ionic  elec¬ 
tronic  conductors  (MIECs)  under  consideration  indicated 
that  most  of  the  ionic  current  was  generated  near  the  cath¬ 
ode/electrolyte  interface  (at  x/L  =  1),  as  shown  in  Fig.  1  for 
a  single  layer  of  LSCuF.  The  approximately  5-8  [xm  region 
of  high  current  generation  is  larger  than  the  usual  three  phase 
boundary  width  discussed  in  the  literature,  probably  because 
of  errors  in  the  estimated  material  parameters  at  this  early 
stage  in  the  analysis.  The  results  are  very  sensitive  to  the 
pore  surface  area/volume  ratio  ( S )  defined  in  Appendix  A. 
Although  the  results  are  consistent  with  LSCuF  performance 
at  low  overpotentials,  a  detectable  amount  of  concentration 
polarization  was  indicated,  which  motivated  the  next  stage 
of  the  analysis. 

2.2.  Material  properties  and  micro  structural  sensitivities 
of  the  model 

The  second  DZA  model  [8]  also  treated  concentration  po¬ 
larization  along  with  the  polarization  due  to  chemical  re¬ 
sistance  (Rch),  and  is  briefly  summarized  in  Appendix  C. 


Fig.  1.  Current  distribution  in  a  homogeneous  LSCuF  cathode. 
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A  second  computer  simulation  code  was  constructed  using 
this  model.  The  code  has  three  main  functions:  (a)  fitting  to 
experimental  V-I  data  to  extract  the  value  of  the  surface  ex¬ 
change  coefficient  directly  from  the  performance  data,  (b) 
variation  of  model  parameters  to  reveal  the  sensitivities  to 
microstructural  variables,  and  (c)  computation  of  the  current 
and  voltage  as  a  function  of  position  in  the  cathode  material. 

Although  the  closed-form  solution  provided  by  this  model 
is  very  useful,  its  limitation  is  that  it  treats  only  homogeneous 
materials.  That  is,  the  microstructural  parameters  may  not 
depend  on  spatial  position.  It  will  be  seen  that  a  piecewise 
numerical  solution  of  the  model  is  required  for  the  final 
design  of  the  cathode. 

Analysis  using  the  extraction  mode  was  applied  to  one  of 
the  best  performing  materials  developed  in  this  laboratory: 
LSCuF.  The  operating  parameters  were  T  =  700  °C  and 
Po2  =  2.12  x  104  Pa  (air)  [1].  The  ionic  and  electronic  con¬ 
ductivities  were  taken  as  0.01  and  lOOS/cm,  respectively. 
The  thickness  was  25  [xm,  porosity  50%,  grain  size  1  [xm, 
pore  size  1  |xm,  with  solid  and  gas  tortuosities  of  1.89  and  2.5 
(see  Appendix  C).  Several  dozen  runs  revealed  that  the  best 
ionic  diffusivity  value  was  5  x  10-8  cm2/s,  which  is  slightly 
higher  than  that  for  LSC  and  LSCoF  (lx  10-8  to  2  x  10~8) 
[6],  but  still  quite  reasonable  in  light  of  the  chemically  ac¬ 
tive  role  copper  plays  in  many  technologies  (e.g.,  getters, 
superconductivity,  etc.).  This  is  probably  due  to  the  elec¬ 
tronic  structure  of  copper,  analysis  of  which  is  outside  the 
scope  of  this  work. 

The  solutions  for  Ks  for  the  various  current  and  voltage 
data  pairs  are  shown  in  Fig.  2.  The  curve  labeled  curs_p  (with 
concentration  polarization)  reproduces  the  data  at  an  ASR 
of  about  0.06  Q  cm2  in  the  ohmic  regime,  and  shows  the  ap¬ 
parent  concentration  polarization  becoming  important  above 
about  0.1  V  as  the  non-polarized  (curs_np)  curve  departs  no¬ 
ticeably  from  curs_p.  When  compared  to  the  solid  material 
(cur_s),  both  curves  show  the  significant  enhancement  pro¬ 
vided  by  the  porous  material.  Lp  (the  characteristic  length  of 
the  three  phase  region)  at  0. 1  V  was  1 .7  |xm  and  Lg  (the  char- 


Fig.  3.  Effect  of  varying  cathode  thickness  with  all  other  parameters  held 
constant:  overvoltage  =  0.1  V. 

acteristic  length  for  concentration  polarization)  was  113  |xm, 
indicating  minimal  concentration  polarization  at  this  volt¬ 
age  (see  Appendix  A  for  definitions  of  parameters).  The  Pq2 
at  the  cathode/electrolyte  interface  had  fallen  to  1.52  Pa  at 
0.1  V.  The  value  of  Ks  at  0.1  V  was  2.17  x  10-6  cm/s,  which 
is  quite  reasonable  compared  to  that  for  LSC  and  LSCoF  at 
700  °C  (3  x  10-6  to  4  x  10-6  cm/s). 

The  above  results  were  then  used  in  the  model’s  sensitivity 
mode  to  study  the  impact  of  varying  microstructural  param¬ 
eters  as  follows.  Fig.  3  shows  the  effect  of  varying  the  cath¬ 
ode  thickness.  It  appears  that  the  present  25  pan  thickness  is 
just  sufficient  to  establish  a  stable  enhancement  of  current 
due  to  the  ion  collection  activity  of  the  outer  regions  of  the 
MIEC  cathode  material.  At  thicknesses  greater  than  25  [xm, 
concentration  polarization  effects  become  noticeable. 

Variations  in  porosity  cause  variations  in  the  pore  surface 
area/volume  ratio  (S)  and  the  tortuosities,  and  are  shown 
in  Fig.  4.  Although  the  enhancement  remains  significant,  it 
is  interesting  that  a  lower  porosity  results  in  a  higher  cur¬ 
rent.  This  may  be  due  to  part  of  the  particular  relation¬ 
ship  (see  Appendix  C)  between  the  ionic  diffusivity  and 
the  surface  exchange  coefficient.  The  characteristic  distance 


Fig.  2.  Analysis  of  LSCuF  data  in  the  extraction  mode.  The  squares  are  the  experimental  data  points.  The  long-dashed  curve  is  the  model  ht  to  the  data. 
Other  symbols  are  explained  in  the  text. 
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Fig.  4.  Effect  of  varying  porosity  at  0.1  V  overvoltage. 


for  the  transition  between  reactivity  and  diffusion  domi¬ 
nance  was  Ld  =  230  pan.  Another  data  set  and  associated 
benchmarking  exercise  may  indeed  give  different  results. 
For  the  present  combination  of  parameters,  at  0.1  V  over¬ 
voltage,  it  appears  that  lower  porosity  may  give  a  slightly 
higher  current.  However,  this  enhancement  is  smaller  than 
that  provided  by  other  microstructural  parameters,  as  shown 
below. 

The  pore  surface  area  decreases  when  grain  diameter  in¬ 
creases,  but  the  tortuosities  are  not  greatly  affected  in  the 
present  formulations.  The  decreasing  pore  surface  area  re¬ 
duces  the  total  ionic  current,  as  shown  in  Fig.  5. 

Fig.  6  shows  that  variations  in  pore  diameter  have  negli¬ 
gible  effect  in  this  regime  because  the  micron- sized  pores 
are  greater  than  the  mean  free  path  of  the  molecules,  thus 
minimizing  the  Knudsen  effect. 

The  conclusion  was  that  smaller  grain  diameters  would 
provide  the  greatest  increase  in  currents,  e.g.,  reducing 
the  grain  size  from  1.0  to  0.25  jxm  could  double  the 
current. 


Fig.  5.  Effect  of  grain  size  on  currents  at  0.1  V  overvoltage. 


pore  diameter,  microns 

Fig.  6.  Effect  of  varying  pore  diameters  at  0.1  V  overvoltage. 


3.  Design  optimization  for  a  two-layer  cathode 

The  above  results  indicated  that  the  most  effective  way 
to  increase  the  current  output  of  the  cathode  was  to  de¬ 
crease  the  grain  size,  while  maintaining  the  porosity  at  a 
high  level.  This  means  that  the  cathode  internal  pore  sur¬ 
face  area  to  cathode  volume  ratio  will  rise  appreciably.  It 
also  means  that  the  cathode  will  be  a  very  fragile  struc¬ 
ture,  and  difficulties  may  be  expected  in  attaching  current 
collectors  to  its  outer  surface.  One  solution  to  this  prob¬ 
lem  is  to  employ  a  two-layer  cathode  system:  a  more  struc¬ 
turally  robust  outer  layer  with  larger  grains  and  an  inner  layer 
(next  to  the  electrolyte)  with  smaller  grains.  The  outer  layer 
should  have  enough  porosity  to  easily  transport  a  steady 
supply  of  air/02  to  the  more  chemically  active  inner  layer. 
The  following  design  is  based  on  the  above  cited  material 
properties  extracted  from  experimental  data  (0.78  A/cm2  at 
0.1  V).  The  only  change  is  that  the  base  porosity  was  re¬ 
duced  to  45%  to  take  advantage  of  the  small  enhancement 
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in  Fig.  4,  and  to  ensure  numerical  stability  in  the  calcula¬ 
tions. 

The  design  of  each  layer  is  described  below.  It  is  impor¬ 
tant  to  note  how  the  model  was  applied  in  these  cases.  Each 
layer  was  treated  as  an  individual  cathode,  i.e.,  the  model  was 
applied  twice,  once  for  the  outer  layer  and  once  for  the  in¬ 
ner  layer.  This  was  necessary  because  of  the  single-material 
limitations  of  the  model.  Thus,  the  results  should  be  con¬ 
sidered  as  approximate  until  they  can  be  confirmed  by  more 
detailed  analysis  with  a  model  that  treats  the  spatial  varia¬ 
tion  of  material  properties  explicitly. 

3.1.  Outer  layer 

The  following  microstructural  parameters  were  de¬ 
fined:  thickness  =  25  p,m,  porosity  45%,  pore  diameter 
=  1  (Jim.  grain  diameter  =  1  pan,  surface/volume  ratio  S  = 
50,000  cm-1,  gas  tortuosity  =  2.22,  solid  tortuosity  =  1.83 
(effective  path  for  ionic  and  electronic  conduction).  Results 
are  shown  in  Fig.  7. 

The  5  |xm  active  region  adjacent  to  the  electrolyte  (at  v  = 
0)  is  evident.  The  Po2  at  the  5  p,m  distance  was  estimated  to 
be  1.62  x  104  Pa,  the  overvoltage  ( [rj )  4.69  x  10~3  V,  and  the 
ionic  current  0.0106  A/cm2.  These  values  define  the  bound¬ 
ary  conditions  for  the  inner  layer,  at  the  inner/outer  layer 
interface.  While  the  overvoltage  and  current  at  the  5  p,m  po¬ 
sition  are  good  numbers,  it  is  important  to  note  that  the  Po2 
value  was  very  approximate.  This  is  because  of  limitations 
in  the  closed-form  solution  of  the  model,  which  did  not  sup¬ 
ply  a  solution  for  the  oxygen  chemical  potential  for  this  par¬ 
ticular  geometry.  It  was  thus  necessary  to  estimate  the  Po2 
using  a  Nernstian  approximation,  which  will  be  improved 
upon  in  later  work. 

Two  designs  for  the  inner  layer  are  discussed  next,  the 
objective  of  which  is  to  increase  the  current  output  above 
that  seen  for  the  single  layer  in  Fig.  7. 


3.2.  Inner  layer 

3.2.1.  Small  grain  diameter  (0.25  /im) 

Because  of  the  smaller  grains,  the  surface/volume  ratio  is 
increased  to  S  =  1.22  x  106cm-1.  The  pore  diameter  and 
gas  tortuosity  were  unchanged,  but  the  solid  tortuosity  is  re¬ 
duced  to  1.67.  Fig.  8  shows  that  it  is  theoretically  possible 
to  increase  the  output  of  the  two-layered  cathode  to  about 
3.1  A/cm2  at  0.078  V.  The  Po2  at  the  cathode/electrolyte  in¬ 
terface  is  about  5.17  x  102Pa  using  the  above  mentioned 
Nernst  estimation  method,  indicating  that  concentration  po¬ 
larization  was  not  dominant.  In  the  figure,  recall  that  curs_np 
is  the  current  with  no  accounting  for  concentration  polariza¬ 
tion,  curs_p  accounts  for  concentration  polarization,  cur_s  is 
the  current  for  a  solid  material,  and  r/  is  the  overvoltage. 

However,  there  is  some  question  that  this  design  may  not 
be  manufacturable:  the  smaller  grains  require  lower  sintering 
temperatures,  and  overlaying  a  larger-grained  20  [xm  thick 
outer  layer  would  require  higher  sintering  temperatures  that 
would  over-sinter  the  inner  layer.  For  this  reason,  a  design 
was  also  generated  for  an  inner  layer  with  a  larger  grain 
size. 

3.2.2.  Larger  grain  diameter  (0.5  \im) 

The  surface/area  ratio  for  this  design  was  62,777  cm-1, 
the  solid  tortuosity  was  1.65,  and  the  Pq2  at  the  cath¬ 
ode/electrolyte  interface  was  less  than  1  x  10-5  Pa  by  the 
Nernstian  approximation.  Because  of  the  low  Po2 ,  the  cal¬ 
culation  actually  failed  at  0.4  jam  from  the  interface.  Fig.  9 
shows  that  the  overvoltage  was  appreciably  higher  (0.12  V), 
although  the  actual  current  output  may  achieve  values  near 
those  of  the  0.25  p,m  grain  case  at  the  interface  proper. 
However,  the  interface  would  probably  be  oxygen-starved, 
thus  eliminating  any  benefit  from  oxygen  impingement  di¬ 
rectly  on  a  chemically  stable  thin  layer  (i.e.,  1  pan  ceria) 
between  the  electrolyte  and  LSCuF  cathode  material. 
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Fig.  7.  Ionic  current  and  overvoltage  for  the  cathode  outer  layer. 
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Fig.  8.  Ionic  current  and  overvoltage  (eta)  for  the  small  grain  inner  layer. 
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Fig.  9.  Ionic  current  and  overvoltage  (eta)  for  the  cathode  inner  layer  with  larger  grain  diameter. 


It  appears  that  the  best  design  is  the  0.25  \xm  grain,  5  p,m 
thick  inner  layer.  However,  the  manufacturability  of  this  de¬ 
sign  must  be  proven  experimentally. 

4.  Conclusions 

The  models  employed  to  this  point  are  suitable  for  a  ho¬ 
mogeneous  material  where  the  properties  and  microstruc¬ 
tures  do  not  vary  spatially.  However,  the  micro  structural 
design  of  the  two-layered  cathode  described  above  has  ex¬ 
tended  these  closed-form  models  appreciably  beyond  their 
intended  range  of  application,  and  should  be  considered  as 
approximate. 

The  candidate  material  LSCuF  investigated  in  this  anal¬ 
ysis  has  exhibited  a  high  affinity  for  oxygen,  and  may  be 
suitable  for  a  cathode  material  if  structural  and  chemical 
stability  can  be  maintained  during  long-term  high  power 


density  operations  at  elevated  operating  temperatures.  Ex¬ 
perimental  efforts  are  underway  to  investigate  the  stabil¬ 
ity  of  this  material,  and  will  be  reported  in  a  separate 
paper. 
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Appendix  A.  Model  parameters 


The  corresponding  current  for  a  porous  material  is  given  by 


The  microstructural  parameters  are: 

•  the  porosity,  0; 

•  the  grain  diameter,  dg\ 

•  the  pore  diameter,  dp; 

•  the  internal  pore  surface  area  per  unit  volume,  S ; 

•  the  thickness  of  the  cathode,  L; 

•  the  solid  tortuosity  for  electronic  or  ionic  conduction,  ts; 

•  the  tortuosity  for  diffusion  of  the  gas  in  the  pores,  tg. 

Operational  parameters  that  can  be  varied: 

•  the  overvoltage,  V; 

•  the  total  gas  pressure  at  the  cathode  surface,  Ptot; 

•  the  Po2  at  the  external  surface  of  the  cathode; 

•  the  operating  temperature,  T. 

Material  properties  input  to  the  model: 

•  the  ionic  diffusivity  for  oxygen,  Aon; 

•  the  electronic  and  ionic  conductivities,  cre  and  o\\ 

•  the  surface  exchange  coefficient,  Ks; 

•  the  ambipolar  diffusion  coefficient,  Ae  =  Aon/(l  + 

•  the  ionic  concentration  in  the  solid  material,  C\\ 

•  the  characteristic  distance  for  the  transition  between  re¬ 
activity  and  diffusion  domination,  Ld  =  Ae/A; 

•  the  characteristic  length  of  the  three  phase  bound¬ 
ary  region,  or  region  of  high  reactivity  near  the  cath¬ 
ode/electrolyte  interface,  Lp; 

•  the  characteristic  length  for  concentration  polarization  ef¬ 
fects,  Lg; 

•  the  chemical  resistance  due  to  ionic  conduction  and  sur¬ 
face  exchange,  Ah- 

Calculated  output  from  the  model: 

•  the  ionic  current  for  a  solid  material,  If ; 

•  the  ionic  current  for  a  porous  material  with  no  concentra¬ 
tion  polarization,  /np; 

•  the  ionic  current  for  a  porous  material  including  concen¬ 
tration  polarization  effects,  /p; 

•  the  ASR  of  the  cathode  material,  given  by  the  current 
divided  by  the  overvoltage. 


Appendix  B 


The  DZA  model  [7]  for  small  concentration  polarization 
is  summarized  as  follows. 

The  ionic  current  for  a  solid  material  is  defined  as 


ACiA/zg 

RT 


(B.l) 


where 


RT 


In 


A)?  external 


A)?  interface  / 


v 


(B.2) 


I\(x)  =  I{MG{x), 


(B.3) 


where  the  material  factor  is 


M  = 


rs 


(B.4) 


and  G  is  a  geometric  factor  given  by 

aexp[(L  —  x)/Lry]  +  exp[v/LD] 

G(x)  = - — - C—Ei - F '  p  ,  (B.5) 

aexp[L/Lp]  +  1 


where  v  is  the  distance  from  the  cathode/electrolyte  inter¬ 
face,  and 

1  -\-  u 


a  = 


1  —u 


(B.6) 


u  = 


1/2 


SLd 


and 


A(  1 


(B.7) 


(B.8) 


Eq.  (B.3)  describes  the  ionic  current  in  the  porous  solid  in 
the  Gerisher  limit,  where  the  oxygen  potential  is  constant. 
The  source  of  this  ionic  current  is  the  oxygen  molecular 
flux  through  the  pores:  the  molecular  current  is  converted  to 
ionic  current  through  the  surface  exchange  mechanism.  The 
molecular  flux  through  the  pores  can  be  expressed  as: 


j  = _ ±  dc§ 

s  rgZ)g  dx 


(B.9) 


where  Dg  is  the  diffusivity  of  the  O2  in  the  pores  (corrected 
for  Knudsen  effects)  and  Cg  is  the  concentration  of  the  O2 
in  the  pores,  which  depends  on  position  v  in  the  porous 
material.  When  all  the  molecular  O2  is  converted  to  ionic 
current,  we  have 

A  —  2^i(max)*  (B.10) 


But  this  does  not  admit  spatial  dependence  of  A)2  • 

To  obtain  that  spatial  dependence,  we  depart  from  the 
original  developments  by  DZA  [7],  as  follows.  Note  that 
at  any  position  x,  the  time  rate  of  change  of  the  molecular 
current  is  proportional  to  the  spatial  gradient  in  the  ionic 
current: 


d/g  _  1  dl\ 

d  t  2  dx 


(B.ll) 


From  the  continuity  equation,  we  also  have 

Pi  =  _d4 

dt  dx 


(B.12) 


Substituting  Eqs.  (B.3),  (B.9)  and  (B.12)  into  Eq.  (B.ll)  re¬ 
sults  in  a  second  order  differential  equation  in  the  oxygen 
partial  pressure  as  a  function  of  spatial  position  x  This  is 
a  two-point  boundary  value  problem,  and  is  integrated  nu¬ 
merically  by  a  relaxation  method.  The  numerical  solution  is 
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Fig.  10.  Benchmark  case  for  LSCoF  at  0.5  V,  L  =25  |xm,  T  =  700  °C. 


then  substituted  into  the  Nernst  equation  (since  concentra¬ 
tion  polarization  is  expected  to  be  negligible)  to  obtain  an 
estimate  of  the  local  overpotential  between  two  points  x(i) 
and  x(i  —  1): 


RT 
4 F 


In 


P02W 0) 

Po2(x(i  -  1)) 


(B.13) 


The  chemical  resistance  Rc h  (or  ASR  for  no  concentration 
polarization)  is  calculated  from 


RT 


2F2K,Q 


ion 


1/2 


(B.14) 


and  the  current  increment  between  x(i)  and  x(i  —  1)  is  cal¬ 
culated  from  Ohm’s  law,  including  the  geometric  factor  G 
for  the  ionic  flux,  by 


A/(*(/))  = 


A  T](x(i)) 

Rch  x  G(x(i)) ' 


(B.15) 


The  current  increments  are  then  added  to  find  the  local 
current  as  a  function  of  position.  The  method  was  bench- 
marked  using  a  case  for  LSCoF  discussed  in  detail  in  Ref. 
[5].  Results  (shown  in  Fig.  10)  indicate  that  the  majority  of 
the  current  is  generated  within  about  2.5  jxm  of  the  cath¬ 
ode/electrolyte  interface  (at  x/L  =  1),  in  reasonably  good 
agreement  with  the  previous  calculation  of  Ld  =  3  pm, 
shown  in  Fig.  8a  of  [5].  Results  for  LSCuF  are  shown  in  the 
main  text. 


Appendix  C 

C.l.  Micro  structural  design  including  concentration 
polarization 

The  second  DZA  model  [8]  solved  two  coupled  differen¬ 
tial  equations  for  the  chemical  potential  /x  of  the  ambipolar 


(ion  +  electron)  species  and  for  the  chemical  potential  /x 
of  the  oxygen  molecules,  as  a  function  of  position  v  in  the 
porous  cathode  material.  The  ionic  current  is  given  by 

,  _  1-0  £>ie  C\  d[i  ^ 

—  r)rr  a  ’ 

rs  RT  ax 

where 


/X  —  /X0 


gll  exp 


+  £12  exp 


+  [£21*  +  £22]^12  -  A/Xg, 


(C.2) 


where  /xo  is  the  ambipolar  chemical  potential  at  the  cath¬ 
ode/electrolyte  interface  (x  =  0),  and  the  gu  are  constants 
determined  from  the  boundary  conditions  The  gu  equa¬ 
tions  are  several  pages  long  and  are  not  repeated  here. 
Several  typographical  errors  in  the  original  publication 
have  been  corrected  for  the  present  work.  The  are 
components  of  the  eigenvector  matrix,  and  Lm  is  given 
by 


i  _  i  1 

TT  ~  Jl  +  77’ 


(C.3) 


where 


=  2  L 


P 


cp  Dg  Cc 


-1 1/2 


_rg  1  —  0  Die  C\  _ 


(C.4) 


is  the  characteristic  length  for  concentration  polarization  to 
occur,  i.e.,  if  this  length  is  smaller  than  the  cathode  thickness 
concentration  polarization  is  important.  A/xg  is  the  molec¬ 
ular  chemical  potential  difference  between  the  two  surfaces 
of  the  cathode. 

It  is  important  to  note  that  mathematical  relationships  be¬ 
tween  the  microstructural  parameters  0,  S ,  rs,  and  rg  have 
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been  derived  for  cubic  systems  in  [10],  and  were  also  em¬ 
ployed  in  this  analysis: 


c  36(1-0) 

*5  —  , 

rg 

(C.5) 

where 

rg  =  \dg, 

(C.6) 

,  1-3  q 

b=  o  2 

i-lo-?) 

(C.7) 

3tt/4[2/9  —  q2(l  —  q/2>)] 

1  0=  ,,  ,3 

(i  -  qy 

(C.8) 

rs  =  -(1-0)(1  -q)  In  (-  — l), 
n  \q  ) 

(C.9) 

1 

J' 

s~  <P' 

(C.10) 

and  q  is  the  fractional  consolidation  of  the  grains,  a  sintering 
parameter. 

This  model  thus  provides  a  closed-form,  simultaneous  so¬ 
lution  for  the  two  polarization  terms  of  interest  in  this  work, 


chemical  and  concentration  polarization.  Although  very  use¬ 
ful  (see  main  text),  its  limitation  is  that  it  treats  only  a  ho¬ 
mogeneous  material.  That  is,  the  microstructural  parameters 
are  not  dependent  on  spatial  position. 
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